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Abstract 

We study the fourth order action of the comoving curvature perturbation in an infla- 
tionary universe in order to understand more systematically the de Sitter limit in nonlinear 
cosmological perturbation theory. We derive the action of the curvature perturbation to 
fourth order in the comoving gauge, and show that it vanishes sufficiently fast in the de Sit- 
ter limit. By studying the de Sitter limit, we then extrapolate to the n'th order action of the 
comoving curvature perturbation and discuss the slow-roll order of the n-point correlation 
function. 
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1 Introduction 



Even though inflation has proven to be a very successful paradigm for the evolution of the 
early universe, we have very limited knowledge of the theory of inflation itself. At present, 
even the energy scale at which inflation takes place is unknown within almost 10 orders of 
magnitude. The cosmological data available has reached a precision, which allow us to test 
the generic predictions from inflation of a flat, adiabatic and gaussian scalar spectrum of 
cosmic microwave background (CMB) anisotropies very well [Tj. On the other hand, the 
data is not yet precise enough to enable us to convincingly find sub-leading deviations from 
the generic predictions, which would enable us to discriminate between different classes of 
inflationary models [2]. 

In order to make a transition to the next level in our understanding, it is important 
to understand more systematically what our theory entails to higher order in perturbation 
theory. A lot of effort has already been put into understanding second order perturbation 
theory [3], non-gaussianities (for a review see e.g.[l]), loop corrections [51 [6j [7J [8j [H [10] , and 
even general nonlinear perturbation theory [TTJ [12]. However, there is one question which 
is of particular interest to us. We are interested in what we can learn about higher order 
perturbation theory by studying the de Sitter limit. Therefore, in order to make a step in 
this direction, in this paper we will study the de Sitter limit of cosmological perturbations 
up to fourth order in two different gauges; the comoving curvature gauge and the uniform 
curvature gauge. 

The action of cosmological perturbations to third order has previously been calculated 
in both gauges in ref. [TB] . and later generalized to include higher dimensional operators 
[TT[ [TB] , and multi field scenarios [TB] . It was found that the third order terms in the action 
are slow-roll suppressed in both gauges and thus vanishes in the de Sitter limit. This is 
expected, since the curvature perturbations should be a pure gauge in the de Sitter limit, 
and no longer be dynamical. One may expect that in the uniform curvature gauge the action 
to all orders higher than two should be slow- roll suppressed in the pure de Sitter limit [16 J, 
but an explicit calculation of the fourth order action in the uniform curvature gauge shows 
that the action of inflaton fluctuations is not slow-roll suppressed at fourth order [6], [TTJ. As 
discussed in [17] . this happens because the gravity is highly nonlinear and the nonlinearity 
in the inflaton fluctuations are controlled by the strength of the gravitational interaction 
rather than the slow-roll parameter. However, when the action is reformulated in terms 
of the comoving curvature perturbation, it becomes slow-roll suppressed. This confirms 
that the curvature perturbation in pure de Sitter space is a pure gauge, which can also be 
understood from the fact that the transformation from inflaton perturbations into comoving 
curvature perturbations becomes singular in the pure de Sitter limit. These are the issues 
we will explore in more explicit details in the present paper. 

In the next section, we calculate explicitly the fourth order action of the comoving curva- 
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ture perturbation. This is the main result of the paper. We show that by a direct application 
of the ADM formalism, one obtains a result which is zeroth order in the slow-roll expan- 
sion. Just like in the third order calculation [13], the zeroth order terms can be eliminated 
by partial integrations to give an action, which is suppressed by one power of the slow-roll 
parameter. However, the leading contribution is proportional to the linear perturbation 
equation, and can be eliminated by a change of variables, leading finally to an action, which 
is suppressed by two powers of the slow-roll parameters. 

By studying the de Sitter limit we extrapolate our findings from fourth order to n'th 
order. We argue that the action of inflaton perturbations in the uniform curvature gauge 
has to be slow-roll suppressed to any odd power, while it will remain unsuppressed to any 
even power in the perturbation expansion. We then review systematically the slow-roll order 
of n-point correlation functions of the comoving curvature perturbations up to n = 4, and 
discuss the extrapolation to any n > 4. As an example, we conjecture that the nonlinearity 
parameters of the 5- and 6-point functions are 0(e 2 ), where e is the first slow-roll parameter. 
The nonlinearity parameters of the 3- and 4-point functions, f^L, t nl, have previously been 
shown to be of order 0(e) [T3"l [T7] . The results for the nonlinearity parameter to n'th order 
is summarized in table 2 at the end of section 3. 

Although our motivation for deriving the action to fourth order in the comoving gauge 
is more conceptual, the result may also have more phenomenological applications. It can 
be used for calculating loop effectl and non-adiabatic enhancements of the tri-spectrum 
beyond the slow-roll approximation |18j . 

The outline of the paper is the following. In section 2 we derive the fourth order action 
of curvature perturbations in the comoving gauge. In section 3 we discuss the de Sitter 
limit and the extrapolation of the results to arbitrary order. In section 4 we summarize and 
conclude on our findings. 

2 The action to fourth order in comoving gauge 

The action of cosmological perturbations during inflation has previously been calculated 
both in the comoving and in the uniform curvature gauge to third order, but only in the 
uniform curvature gauge to fourth order. The two gauges have different benefits, and it is 
useful to complete the picture to fourth order. Thus, we proceed by computing the action 
of curvature perturbations to fourth order in the comoving gauge. 

In order to calculate the action, it is convenient to use the ADM formalism [21], and 
slice the 4-geometry into a sandwich structure along the time-coordinate. One can produce 
a series of spatial 3-geometries of uniform time, separated by time steps dt. To connect two 

*In ref. [5] the fourth order action in the uniform curvature gauge was used to demonstrate that the 
loop corrections to inflation can be large, and may have observationally important consequences for the 
predicted tensor-to-scalar relation in models of chaotic inflation. This is consistent with ref. [3 [7] , and has 
later been verified using a combination of the fourth order action in the uniform curvature gauge and the 
8 iV-approach [9l [10] . It is possible to make use of the <5 TV-approach, because the IR contributions to the 
loops also can be described by a particular classical approximation (see also [20]). 
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spatial slices, one introduces the lapse function, N. Similarly one can write the difference 
between the position in the two 3-geometries in terms of the shift-vector N l [21] 



<ew = initial ~ N*dt . (1) 

The quantities A" and N l will later be expanded in series of the perturbation parameters. 
With these definitions, one can write the line element on the ADM form [21], in terms of N, 
N\ 

ds 2 = -N 2 dt 2 + hij (dx i + N l dt) (dx j + N j dt) , (2) 

with being a purely spatial metric. Using this, the Einstein-Hilbert action for a single 
scalar field 
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S = I d A xy/^g 

can be written as 



R+-d^dy-V((j)) 



(3) 



S = \J dtd 3 xVh NK 3 - 2NV{(j>) + N~ l (EijE ij - E 2 ) + N-\<p - N l d v 
- Nh ij 'difidjcp 



(4) 



with IZ3 being the curvature scalar associated with hij. The tensor Ey, which is closely 
related to the extrinsic curvature, is defined as a linear combination of the time derivative 



of hij and the covariant derivative of the shift vector 



E i j = ^{h i j-V i Nj-VjN i ) , (5) 



and E = E\ is the trace of the tensor Eij. 



2.1 Gauge choice, constraints, and solutions 

In this and the following subsection, we will work in the comoving curvature gauge, defined 



as 



= <t>c{t) hj = a{tf{e 2 ^5ij + 7<i ) , 7 « = d aij = 0. (6) 

For the purpose of this paper, we can focus on scalar perturbations of the metric, thus letting 
jij being identical to 0. 

In the present gauge the curvature scalar becomes IZ3 = —4d 2 ( — 2 (<9C)tl an d from the 
action in eq. (j3J), one can derive the equations of motion A^, 

VilN-^E'j-SiE)] = 0, (7) 

and for A^ 

TZ 3 -2V - X 2 (/•-'-, /•-'" - E 2 ) - N- 2 <j> 2 = . (8) 



t"We use the notation d 2 C = h^didjC and (d() = h t; >di(dj(. 
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If one furthermore expands the scalar and shift parameter in terms of three perturbation 
parameters to an arbitrary order, one obtains 

N = 1 + + « (2) + . . . ) N t = d t ( X w + X (2) + ..-) + (A' (1) + A (2) + • • • ) (9) 

with the constraint (cf. Helmholtz theorem) 

&fc = . (10) 

It is now possible to solve equations ([7]) and © order by order for all three perturbation 
parameters. As it will turn out, we only need to calculate the first and second order per- 
turbation parameters. The fourth order terms only appears in the action multiplied by the 
equation of motion for the background field, C , while the third order terms cancel one an- 
other. From eq. ([7]), by taking the divergence and employing eq. ffTOl) . one finds to first order 
in a [15] . 



a 



H-\ • (11) 

Inserting this back into the first order contribution to eq. (J7J, and solving for the first order 
vector perturbation, one also obtains [33] 

= o . (12) 

In eq. (1111) we have introduced the parameter H = a/a. 

Utilizing the same techniques for the second order contributions, we get 



a (2) = ^d~ 2 {& (^a (1) 9 2 x (1) - d i a w d i d j x ll) ) ~ 28 i d j (d j x (1) dX) 
+ & {dXd'd^) + & {d£d 2 x {1) ) 
for the scalar part, and 

/V 2) = 2d~ 2 \d,ja m d 2 X w - dtoP&djX™ + dj(d 2 X w - 2Hd ja {2) 
- di (d jX w dX + d l x w d 3 C) + dXd l d jX (1) 



(13) 



(14) 



for the vector perturbation. In the equations above is the first order part of the third 
perturbation parameter, which we will solve for in a moment. From this we see, that a to 
the lowest order only describe how the metric changes in time, as could be expected from 
its role in the metric. 

Similarly we can solve eq. (|HJ) for the x parameters, whereby one gets 

x (1) = -| + e , 9 2 e = ^ 2 c (is) 
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to the first order [13]. Likewise to the second order in x, we obtain (for readability we refrain 
from inserting the values of the a and f3) 



4Hd 2 X {2) = -2 (d() 2 + 2(<p 2 - 6H 2 )a {2) + 4(d 2 X (1) + 40 2 a (1) C 

- AHd lX w d l C - 3a (1)2 2 + d 2 X w d 2 x w - ^o^ 1 ' ' ° ' 

It is seen, that the j3 and x parameters of a given order depends on the a parameter of 
the same order, as opposed to the a that only depends on parameters of lower orders. This 
translates into a dependence of the lapse function in the shift parameter, which is reasonable 
as the lapse function is linked to the thickness of the time-slices, and in that way affect the 
size of the shift. 



2.2 Fourth order action 

Truncating the action at the fourth order and exploiting the equation of motion for the 
background field 

2 - 6H 2 + 2V = , (17) 

we write the action in terms of the perturbation parameters, and simplify it by using partial 
spatial integrations. This gives 



s w 



dtd'xa 3 { - -Cd 2 ( - 2a m ((d i (d l ( + ( 2 d 2 + <f) 2 c a 



V - 3C« (1) + a (1)2 



~C + C« (1) + a 



„(D 2 



[didjx^d^x^ ~ d 2 X m d 2 x m ] + {QH 2 - 2 )« 



2^,(2) 2 



2[C + a w ] [d i d j x w &d j X™ ~ d 2 X m d 2 x (2) - 2d i d ]X w d i x w d j (] 
2 [2d t d jX {2) d l X w d>( + 2d i d jX m d i X m d j C - d 3 x m d l Cd l X w & (\ 
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;(2) 



(18) 



Though fairly compact in its presentation, it is not clear from this form that the action 
in eq. ([TBI does indeed vanish in the slow roll limit. The zeroth order slow-roll terms can 
however be removed by partial integrations. To make this more evident, it is simpler to 
proceed with finding a gauge transformation between the uniform curvature gauge and the 
present gauge, in order to calculate the action by a gauge transformation of the equivalent 
action in the uniform curvature gauge. 



2.3 Gauge transformation 

One can find the gauge transformation by doing a translation in the time coordinate in order 
to transform between the uniform curvature gauge and the comoving gauge. Below we follow 
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a procedure, which is similar to the one in ref. [13]. The uniform curvature gauge is given 
by 

0(t, x) = <j) c (t) + 5<j)(t, x) , h i:j = aitfdij , (19) 

where we have defined the background field C = (0 |0| 0), such that the tadpole condition, 
(0 \5<j>\ 0) = 0, is satisfied. 

Choosing the time translation given by the vector £ M = (T, 0,0,0), we can write the 
transformed field fluctuation [22] 

(20) 



n! 

n=l 



The sum is nothing more than the Taylor expansion along a vector. Requirering the time 
translation to bring the coordinates to the comoving gauge, we fix £ M by <50(a;^ + £ M ) = 0. 
Evolving order by order we then get for £ = T (1) + T (2) + . . . . 

T« = -** (21) 



m = ^ _ 1 0c 2 (22) 

= ±tl80<l> ~ - ^ - {-^- ~ W (23) 

204 2 03 ^ 3 ^2 05 6 0^ ^ 

As a method for finding the relation between £ in the comoving gauge and 50 in the uniform 
curvature gauge, we perform coordinate transformation above in the metric associated with 
the uniform curvature gauge (gffi) and equate it with metric in the comoving gauge {g^v^)- 
The spatial part of the metric gives after the time translation (again ignoring tensor 
contributions) 

/ig m) = a 2 (t + T)8ij - N^diTdjT + diTNj + djTN, , (24) 

with iV and Ni defined as in eq. (Q. Solving the first order contribution in the above equation 
gives ( = HT or [13] 

C = x = C« ■ (25) 

0c 

No spatial reparametrization is needed for the first order. To higher orders one must do a 
transformation X{ — > Xi + Vi, such that 



- N 2 diTdjT + (diTNj + djTN,) + d iVj + d jVi = exp ^2 ^ t n j 



a 2 (t)^- (26) 



with L n being a parameter of the n'th order. Recalling that L\ = and writing a(t) = e p ^ 
like in ref. |13j . we can write the equation for £ up to third order 

£ = p(t + T)- p(t) + (l 2 + l 3 ) (27) 
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By taking the trace and d l d^ of eq. (I26I) . one can solve to second order 

4i 2 = 2d i T (1) d lX<t> w - 2d- 2 d i d j {d i T {1) d j X4, m ) ~ (<9 4 T (1) <9 4 T (1) - d^d^^d-T^d^)) (28) 
as = 0. Thereby giving [13] 



C = C - / 2 (Cn) =C + \-i^C + \§- 2 C + ^ + 5#C»W - *3-'a i d>(a i c n d j xS») 



1 

— I 

2 



4F 2 



(diCn&Cn-drWdHdiCndjC*)) 



(29) 



where we implicitly defined /^(Cn)- We refer to appendix [A] for the definitions of the functions 
of the type , an d ■ Computing the third order terms in a similar fashion, we get 



2. 3 =d i T^da/ a) + d l T m (3^ 2) + d i T w diX^ ll) - 8r % &&{d i T™d j x^) 

- d~ 2 & {d 1 djT {1) fl^) - d^d'd^diT^djX^) 

- [<9,T (1) <9 4 T (2) - d~ 2 d i d i {d i T^d j T^)\ 

- a (1) [^T (1) ^T (1) - drWdHdiTVdjT™)] 

Repeating the calculation leading to eq. ( |29|) . we arrive at 

C =Cn - h(Cn) ~ h(Cn) 
=(n — MCti) 

6 # 3 ^ n 3 c # 2 ^ 4# 4 ^ 2j) c H n ^ H 3 ^ n H 2 



(30) 



<Pc . „_i i ? c > 



2H 2 



+ 



Cn 



(31) 



2# 2 \ 2# 2 

d- 2 d i d j 



^-dx n d t d-\ n - H-'dXndiCn ) - —d~ 2 && {d£ n d m 



,(2)> 



2 



^ 2 Cn + Cn + Cn 



2# 2 



diCndjd 2 (n - H 1 d i ( n dj( n 



Cn 



2H 2 \ 2H 2 



OjL, „<),<) Cn - H diCndjCn 



W d- 2 d^{d%Cn^) 



+ iff [^Cn - a- 2 a^(^Cn^Cn)] , 
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where we have implicitly defined the third order contribution in terms of fsiCn)- Again we 
refer to appendix [A] for the definition of functions of the type ct§\ x\ , > an d Pij, ■ 

By inverting eq. (13T1) to third order and recalling the definition of Q n from eq. (125]) . we are 
able to calculate the fourth order action of the metric perturbations in the comoving gauge 
from the action of the inflaton field fluctuations in the uniform curvature gauge. Since the 
calculation is trivial but tedious, we will not repeat it here in all details. However, from the 
calculation, one would obtain an action for ( which is suppressed by only one power of the 
slow-roll parameter. This can be seen from the gauge transformation above, which contains 
terms in fs of order but with no slow-roll suppression. To third order in the gauge 
transformation we can just replace ( n with ( in ^3, and applied in the second order action 
for 8(f) we see immediately that the fourth order action in ( will receive contributions, which 
are only suppressed by one slow-roll order. From this one could easily be led to believe, that 
the trispectrum of the metric perturbation ( is of the wrong order in the slow-roll parameter 
e. 

However, since all the leading order slow-roll terms comes from using / 3 in the second 
order action in 8(f), by partial integrations they can all be rewritten as terms in the action 
proportional to the linear equation of motion of perturbations. This complication can be 
avoided by choosing the gauge of ( n , then calculate and afterwards obtaining the 

super- horizon value of (£ 4 ) by virtue of the super- horizon limit of eq. (I3T]) . analogous to the 
calculation in ref. [13]. Once the terms proportional to the linear perturbation equation of 
motion have been eliminated by a change of variables ( — > ( n , one obtains an action in ( n , 
which is suppressed by two powers of the slow-roll parameter. The full action in the new 
variable ( n has the following form 



6 an 2 ■' 
- 2e H -Hr] H ( n + (, 

+ aJ 2) \-6H 2 + & 



c V 2 



( ( 2 ) 



v {^jV" + 2#V4 - H 2 n H e H + 3H 2 e H (r] H - e H ))Q 



(32) 



d^Wd^Cn + 2e H (2(e H - r) H )H( n + ^m„<)'<) \ 



where we have introduced the slow roll parameters associated with the Hamilton- Jacobi 
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formalism: 



e H = 2M/(^g] 2 =— 2 tr, (33) 
u \ //,,„ ) 2M p i 2 H 2 v 1 



2H"{4>) 



1H= 2M pl ^ = -^. (34) 



We have further more added a second order parameter 



2 H"'(4,)H'(4>) , 2 



To leading order they are related to the usual slow-roll parameters by e# ~ e, tjh ~ t) — e. 



3 The de Sitter limit 

Since the curvature perturbation is a pure gauge mode in de Sitter space, the action of the 
comoving curvature perturbation ( must be slow-roll suppressed to all orders in perturbations 
and vanish in the de Sitter limit. However, this does not explain why the third and fourth 
order action of ( n are actually suppressed by the same power of the slow-roll parameter. If 
we want to estimate the slow-roll order of the n-point function of (, we would overestimate 
it, if we only assumed that the action of ( has to be slow-roll suppressed. 

In fact, the slow-roll order of the action is more easily understood in the uniform curvature 
gauge. In this gauge, it has been shown that the third order action of the inflaton field 
fluctuations, S^dcj)), is suppressed by the square root of the slow- roll parameters [15] . while 
the second, and fourth order actions, 5*2(50), S^(5(p), are unsuppressed [6l [17]. Thus, we 
may wonder what requires the third order action, S3 (50), to be slow- roll suppressed, when 
there is nothing which forces the second and fourth order actions to be slow-roll suppressed. 

However, if there are third order terms of 5(f) in the action, which are unsuppressed and 
survive in the pure de Sitter limit, they would indicate an instability of the classical de 
Sitter vacuum. As an example, we will study a toy unsuppressed third order term. In the 
fourth order action there are unsuppressed terms of the type 5(f> d~ 2 (did <j)d l 5(f)) etc. Let us 
assume for a moment that to third order we have similar unsuppressed terms induced by 
perturbations of the metric. Then the interaction Hamiltonian for the perturbations will 
take the form 

H I = J d 3 ya 3 5(p(4) c + 3H<j) c + V')+gO(5(f) 3 ) + ... , (36) 

where the unsuppressed O(<50 3 ) toy term could be any operator of the type 50 3 , 8(f) S(f>, 
5(f>d~ 2 (di5(f)d' l 5(f)), etc. The dots represents any terms to higher order in slow-roll and per- 
turbation theory, and g is a coupling constant. 

The term 0(5(f> 3 ) will give a contribution to the tadpole diagram in fig. (1), which will 
lead to a one-loop correction to the equation of motion of the classical background field. The 
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tadpole condition yields 



= (5(f)} =4> c + 3H<j> c + V' + gT t , (37) 

where gT t denotes the amputated tadpole contribution, and gives the one-loop correction to 
the background equation of motion. In the simplest case of a massless scalar field and with 
0(5(p 3 ) = H5<f) 3 , the tadpole contribution would become 

gT t = 3gH(6<j ) 2 ) = ^- 2 HH. (38) 

In this case, the time-independent de Sitter solution is destabilized by the tadpole. If the 
toy term 0(5</> 3 ) has a more complicated form involving derivatives, the infrared divergency 
will be absent^, and the tadpole contribution will not grow indefinitely, but rather approach 
a constant. However, in the case of a massless scaler field in de Sitter, a time-independent 
tadpole contribution gT t = const, will give an effective linear contribution to the potential of 
the background term, similar to a source term, which will yield the potential of the massless 
field unbounded from below. In fact, the solution to eq. (|37|) . with V = and gY t = const. 
at late times is 

= -§^, (39) 

which is inconsistent with a time-independent de Sitter solution. Thus, if unsuppressed 
third order terms were allowed, classical de Sitter space with a massless scalar field would 
be destabilized. 




Figure 1: Tadpole diagram 



Since any unsuppressed odd order terms in the action would lead to a non- vanishing 
tadpole contribution, we conclude that all odd order terms in the action should be slow-roll 
suppressed, while even order terms are not slow-roll suppressed. We can then extrapolate 
our results from n < 4 to any n, as shown in table 1. Using the extrapolation in table 1, 
we can calculate the order of magnitude of any n-point correlation function in single field 
inflation, as shown in table 2. We can also generalize the nonlinearity parameter /jvl to n'th 
order up to a numerical factor of order one 

(O « Mt 1 > ( 4 °) 

■^The actual tadpole was calculated to leading order in slow-roll in ref. [10]. Note, that to leading order 
in slow-roll, the IR divergent terms are not present. The IR divergent contributions to the tadpole, which 
will dominate the tadpole contribution at late times, appears to higher order in slow-roll, and can be found 
by taking the appropriate infrared limit in the action as in ref. [H [8] . 
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where V( is the power-spectrum of comoving curvature perturbations. To second order 

(3) 

the generalized nonlinearity parameter coincides with the usual one f NL = Jnl for the bi- 
spectrum, which was initially calculated in ref. (3j [13]. To third order it coincides with the 
nonlinearity parameter for the tri-spectrum fjp L = tjvl, which was calculated in ref. [T7] . 



Order 


S(8<f>) 


S(Q 


S(C) 


2nd 


0(1) 


0(6) 


0(e) 


3rd 


0(e 1 ' 2 ) 


0(e 2 ) 


0(e) 


4rd 


0(1) 


0(e 2 ) 


0(e) 



2nth 



[ j2n"£l)th |t "0(e^J TO(e I 0(e)" ! 
Table 1: Slow- roll order of the action to n'th order. 



0(1) 



p 


(W 


<C P > 




2 


0(H 2 ) 


0(e- 1 if 2 ) 


0(1) 


3 


0{e l / 2 H A ) 


0(e- 1 ^ 4 ) 


0(e) 


4 


0(H 6 ) 


0(e- 2 iJ 6 ) 


0(e) 


2n 


0(H 2p ~ 2 ) 


0(e' p/2 H 2p ~ 2 ) 


0(e p / 2 " 1 ) 


2n + 1 


\0{e^ 2 H 2 v- 2 ) 


0{e^-^ 2 H 2 P- 2 )) 


o( e r P -ir/2)- 



Table 2: Slow-roll order of the n-point functions and generalized nonlinearity parameter. 

It is easy to verify that the action S(5<p) to any even order, 2n, in perturbations will be 
unsuppressed in the slow-roll parameters, as it will contain contributions from a (2)2n , which 
is unsuppressed in the slow-roll parameters. Similarly, to odd orders a {2) will always appear 
in combination with some a {n) or x (n) to odd order, say a w which is slow-roll suppressed. 

As an example we predict that the nonlinearity parameter related to the 5- and 6-point 
function is = f£ ] L = e 2 . 

Finally, let us briefly discuss what would be the effect of including gravitational wave 
modes in the analysis. For gravitational waves the discussion of the de Sitter limit is a little 
more involved. One can have unsuppressed odd n-terms in the action with gravitational 
waves, jij, of the form diS^djScjyyV , since they will not contribute to the tadpole of the scalar 
field fluctuation. However, terms like Scjrfijj 1 -' have to be slow-roll suppressed, because they 
will contribute to the tadpole with a graviton circulating in the loop. This agrees with 
the results of ref. [T3], where the third order action including gravitational waves has been 
calculated. 
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4 Conclusions 



We have calculated the fourth order action of the comoving curvature perturbation, and the 
third order gauge transformation between the uniform curvature gauge and the comoving 
gauge. We have shown that the fourth order action of the comoving curvature perturbation 
is suppressed by one power of the slow-roll parameter, and when terms proportional to 
the linear perturbation equation is removed by a field redefinition, the relevant action is 
suppressed by two orders in slow-roll. This is consistent with the previous calculations in 
the uniform curvature gauge. 

In the uniform curvature gauge, it has been shown that the action of inflaton field per- 
turbations is unsuppressed to fourth order, even though the action is slow-roll suppressed to 
third order. We argued in the previous section, that to any odd order the action of inflaton 
field perturbations has to be slow-roll suppressed in order for the classical tree-level de Sit- 
ter vacuum to be stable, since otherwise a non-vanishing tadpole contributions would shift 
the de Sitter solution. Since the even order terms in the action does not contribute to the 
tadpole, they do not have any destabilizing effect on the background, and are thus allowed 
to survive the de Sitter limit. 

Extrapolating this argument to n'th order, we can estimate the slow-roll order of the 
action to any order, and use it to estimate the n-point function of inflaton perturbations to 
arbitrary order. Although the non-gaussianity from the bi-spectrum is difficult to measure 
at present, and the higher order sources of non-gaussianity will be even harder to detect, 
the result gives a useful theoretical insight that can be used for different purposes. If one is 
interested in possible large non-gaussianities from more exotic models of inflation, it is useful 
to understand carefully the simplest case of single field slow- roll inflation first. In fact the 
calculated fourth order action of comoving curvature perturbations can be used to calculate 
the enhancement of the tri-spectrum due to non-adiabatic effects during inflation [18] . 
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